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Abstract 

> 

' Various processes that are forbidden in the vacuum due to angular momentum 

conservation can occur in a medium that is isotropic and does not carry any angular 
momentum. We illustrate this by considering explicitly two examples. The first one 
is the decay of a spin-0 particle into a photon and another spin-0 particle, using a 
model involving the Yukawa interactions of the scalar particles with a charged fcrmion 
field. The second one involves the decay of a neutrino into another neutrino and a 
graviton, in the standard model of particle interactions augmented with the linearized 
gravitational couplings. 



1 Introduction 

Some physical processes that are not allowed to occur in the vacuum can occur in the 
presence of a background medium. We can classify such processes into two broad classes. 
The first one refers to those processes that are disallowed in the vacuum for kinematic 
reasons. One classic example of this type is the Cerenkov radiation of a charged particle, 
which can occur in a medium due to the fact that the photon dispersion relation is modified 
by the background effects. We are not concerned with processes of this type here. 

The second class consists of processes that do not occur because the transition matrix 
element is zero in the vacuum. In general, whenever that is the case, it can be attributed 
to some conservation laws, which are in turn the consequences of the symmetries of the 
Lagrangian. It is common to refer to such process as being forbidden. Thus, when the 
medium is not invariant under the full symmetry group of the Lagrangian, the relevant 
transition matrix element can be non-zero when the background effects are taken into 
account. 
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In a recent paper [T] we considered a subset of processes in this class, namely those 
which are forbidden by helicity arguments, or angular momentum conservation. We specif- 
ically considered the radiative decay of a spin-0 particle into another spin-0 particle, the 
decay of a spin-1 particle into two photons, the gravitational decay of a spin-0 particle 
into another spin-0 particle and the gravitational decay of a spin- 1/2 particle into another 
spin-1/2. By performing a form-factor analysis in each case, we reviewed the arguments 
that show that the amplitude for the process in the vacuum vanishes, and then demon- 
strated that the amplitude need not vanish if the process occurs in a medium, even if 
the medium is homogeneous and isotropic and therefore does not carry any net angular 
momentum. 

The aim of the present work is to pursue this further to confirm that this is indeed 
the case, by computing the amplitude for such processes in viable models and verifying 
that it does not vanish due to some unexpected reason. Here we consider the radiative 
decay of a scalar particle and the gravitational decay of a neutrino as illustrative examples 
of that kind of process. In the scalar decay case, the model consists of two electrically 
neutral scalar fields coupled to a charged fermion field via Yukawa interactions, and the 
medium is assumed to consist of a thermal background of the charged fermions. In the 
neutrino case, the medium is a thermal background of electrons and, in order to include the 
graviton interactions, the Standard Model couplings are supplemented with the linearized 
gravitational couplings of the particles involved. In either case, the thermal background 
is parametrized by the Fermi-Dirac momentum distribution functions of the background 
particles in the usual way. 

The fact that these processes are forbidden in the vacuum by angular momentum 
conservation arguments, implies that, in the medium, their angular distribution and dif- 
ferential decay rates have a distinctive form. This could lead to observable consequences in 
specific physical contexts despite the fact that there may exist other competing processes 
with comparable total rates. 

In general, the presence of the medium modifies the dispersion relations and wave 
function normalization of the particles that participate in the process. Those corrections 
affect the kinematics, but their relative importance depends on the particular application 
and the physical context of the calculation. Here we assume that the situation is such 
that those corrections are negligible. However, the inclusion of those corrections in the 
calculation of the rates is straightforward [2j, and it should be kept in mind that they may 
be important and need to be included in specific applications. 

In Sec. [2] we consider in detail the scalar radiative decay. There we summarize the the 
form-factor analysis presented in Ref. pQ, show that the on-shell amplitude depends only 
on one form factor, and the expressions for the total and the differential decay rates in 
terms of the form factor are given. We then carry out the calculation of the amplitude 
in a simple model involving a background of charged fermions, and the one-loop formula 
for the on-shell form factor is obtained in terms of integrals over the background fermion 
distribution functions. The integrals are evaluated explicitly for some particular cases of 
the distribution functions. The analogous calculations for the flavor-changing gravitational 
decay of a neutrino are presented in Sec.[3j Seed] contains some general and concluding 
remarks. 
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2 Radiative decay of a spinless particle 



2.1 Kinematical considerations 

In this section, we consider the process 

0(p)- W)+7(<r), (2-1) 

where, <f> and <f>' denote the scalar (spin-0) particles, 7 denotes the photon, and p, p' and 
q denote the corresponding momentum vectors. As mentioned in the Introduction, we 
neglect the effects of the medium on the dispersion relations and wave function normaliza- 
tions in the calculation of the decay rate. While it should be kept in mind that in general 
those corrections must be taken into account in specific applications, they are not essential 
for our purposes here. Moreover, they can be included in the calculations that follow in a 
straightforward way if needed. Therefore for our purposes, we assume the vacuum on-shell 
relations 

2 2 
p = m , 

p' 2 = m' 2 , 

q 2 = 0. (2.2) 

The medium is assumed to consist of a homogenous and isotropic thermal background 
of particles. Besides the thermodynamic variables such as temperature and chemical 
potentials, such a medium is characterized by its velocity four vector v^j3iHl- 
The amplitude can be written in the form 

M = e*»(q)j ll , (2.3) 

where e^(q) is the photon polarization vector which satisfies 

<?%(<?) = 0, (2.4) 

and jn is the matrix element of the electromagnetic current, which satisfies the transver- 
sality condition 

^ = 0. (2.5) 

In general is a function of p^ and q^, which are the only independent momenta in 
the problem. When the process takes place in a medium, can also depend on v^. In 
Ref. [T], it was pointed out that the most general form for the on-shell vertex function 
subject to the transversality condition of Eq. (|2.5p is 

j IM = a\p-qv'*-q-v P >*] + 6e^ 7 p a A 7 • (2.6) 

This implies that the decay amplitude can be written as 

M = (aF; u + bF; u y p \ (2.7) 
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where 



Fftv = tfxQv ~ Q^u , (2.8) 

and = ^s^apF 0113 is the dual. The form factors a and b appearing in the amplitude 
are Lorentz invariant functions of p, q and v. 

The term associated with the form factor b is parity violating due to the presence of 
the Levi-Civita tensor. In a system in which the interactions are parity conserving and 
the background medium is parity symmetric, i.e., is not spin-polarized, the form factor b 
is zero.. In the calculation that we carry out subsequently in this paper, we assume that 
this is the case, and therefore we set b = and write 

M = aF; u v»p v . (2.9) 

Without loss of generality, we will take the form factor a to be real. 

Choosing the convention in which the z-axis points in the direction of the motion of 
the decaying particle, the components of the vectors and p^ in the rest frame of the 
medium are 

^ = (1,0), pP = (E,Pz), (2.10) 

and we will denote by 

V = P/E (2.11) 
the velocity of the decaying particle. The differential decay rate is then 

dT u (1 - F 2 ) 3 / 2 

d{cos9) ~ 16vrm 2 (1 - VcosO) 2 

where 6 is the angle between q and the z-axis, 



\M\\ (2.12) 



2 II 

m — m , 
"0 = ^-, (2.13) 

and 

W = ElM| 2 > (2.14) 

pol 

with the sum being over the two polarization states of the photon. It may be convenient 
to express the differential rate in terms of the photon energy, which is given by 

jYi 2 — ?7l' 2 

" = 2£(i-ycos*r (2 - 15) 

Thus, 

dT _ 1 (1 - V 2 



dto 16irm 2 V 



■\M\ 2 , (2.16) 
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and the total rate is 



167T771 2 V 

where 



r 77 j / <M-M| 2 , (2-17) 



1 + V ' 



(2.18) 



and wo has been defined in Eq. (|2,13l) . 

For the present case, in which the amplitude is given by Eq. (|2.9 



pol 

= a 2 P 2 u? sin 2 9 , (2.19) 

where we should remember that uj and 9 are related by Eq. (|2.15|) . Thus, using Eqs. (|2.12|) . 
(I2TT3]) and (|2TT5l) 

dT a 2 ( m 2 — m' 2 \ rr n, Tr 9 N q /9 sin 2 9 

V 2 (l-V 2 ) 3/2 - —j. (2.20) 



d(cos6») 16ir\ 2m j v ' (l-l/cosfl) 4 ' 

The integration that remains to obtain the total rate, either from Eq. (|2.20p or Eq. 
(|2.17|) . cannot be performed until we have the explicit formulas for the form factor, since 
in general it depends on 9 or, equivalently, u. Here we observe that, since the functional 
dependence of the form factor on 9 will be different depending on conditions of the fermion 
background, so will be the angular distribution. 



2.2 The model and the diagrams 

In order to perform calculations to evaluate the form factors we consider a model containing 
two neutral scalar fields 4> and eft', and a charged fermion /, with the interaction Lagrangian 

= -A/70 - X'fM , (2-21) 

in addition to the standard electromagnetic coupling of /. The medium is assumed to be 
a thermal background of the fermions /. The lowest order diagrams are shown in Fig. [H 
where the internal line in the loops represent the fermions in the thermal background. We 
write the fermion thermal propagator in the form 

iS(l) = %Qf+mf)T(J), (2.22) 

where 

T(l) = A(0 - 2ni5(l 2 - m 2 f )ri{l) , (2.23) 
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l + q 




(a) 




Figure 1: One-loop diagrams for the process of Eq. (|2.ip in a thermal background of 
fermions. Each line is labeled with its corresponding momentum variable. 



with 



A(0 



I 2 — m 2 f + it 



(2.24) 



and 



77(0 



Q(l-v) 



+ 



e(-z-u) 



P(lv-n f ) _|_ ^ e -0(l-v-Hf) _|_ i ' 



(2.25) 



/? and fJLf being the inverse temperature and chemical potential of the fermion gas, respec- 
tively. 

The contributions to from the two diagrams are given by 



i4e f XX' 



dH 



LJl + q,l+ p, l)T(l + q)T(l + p)T(l), 



(2tt) 4 ^ 

He/XX' I j^L^ I -P,l- q)T(l)T(l - p)T(l - q) 



where ef is the electric charge of the fermion /, and 



(2.26) 
(2.27) 



Lp(pi,P2,P3) = jT^hii(J>i + mf)(j>2 + m f )(j> 3 + m f )] , 

= Pl^P2 • P3 ~ P2^PZ ■ Pi + PZtlPl ■ P2 + {Pin + P2 M + P3n) m } • ( 2 - 28 ) 

This expression confirms the statement made in Sec. 12. 1\ that the form factor b that 
appears in Eq. (|2.6p vanishes in this model. As stated there, this is a consequence of 
parity invariance of the interactions of Eq. (|2.21j) and the polarization independence of 
the fermion distribution functions in the medium. 
Using the relationships 



LJ-Pl,-P2,-P3) 



L^(p3,P2,Pl) , 



(2.29) 
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Eqs. (|2,26p and (|2,27p can be combined in the form 



f dH 

jfi = i4e f \\ J ——L^l + qJ+pJ) 

x [T(l + q)T(l + p)T(l) - T(-l - q)T(-l - p)T(-l)] . (2.30) 

When Eq. (|2.22p is substituted into Eq. (|2.30p various terms are produced. The purely 
vacuum terms, which do not contain any factor involving rj(l), give zero since A(— I) = 
A (7). The terms containing three factors of r\ also yield zero due to the on-shell conditions 
implied by the delta function in Eq. (|2.25p . The terms containing two factors of rj will 
give an absorptive contribution to the amplitude. Here we will assume that the kinematic 
regime is such that those absorptive terms are zero, i.e., that the initial (f> state is below 
the // pair production threshold. Therefore the only the terms that survive are those 
that contain one factor of 77, and by making appropriate shifts of the integration variables 
in some of those terms, we obtain 

U = 4 e/ AA' J ^Ls(l 2 -m})[7 1 (l)-7 1 (-l)][L^l + q,l+p,l)A(l+p)A(l + q) 



(27T) 

+L^l -p + q,l,l- p)A{l - p)A{l -p + 
+L fl (l, l+p-q,l- q)A{l - q)A{l +p- 



(2.31) 



2.3 Transversality condition 

The transversality condition q^ = is easily verified explicitly with the help of the 
identity 

(Pi -P3) f *L fl (p l ,p 2 ,p 3 ) = { P2 -p3 + m})A- 1 (p 1 ) - (p 2 -pi + m})A- 1 (p 3 ), (2.32) 

which follows simply from Eq. (|2.28|) . Using this identity to rewrite the various factors of 
q^L^ that appear when we contract Eq. (|2.3ip with q^, we obtain 



h = 4e\\' f -^8(l 2 -m})[ V (l)- V (-l)}{A(l + p)ll.(l+p)+m 2 f ] 



(2vr)= 

+ A(l -p)[l-{l-p)+ mj] - A(l-p + q)[l-(l-p + q) + 



m 



- A(l+p-q)[l-(l+p-q)+m}}} , (2.33) 

where we have also used the relation 

5(l 2 -m 2 f )A- 1 (l 2 -m}) = 0. (2.34) 

Since the factor 7](l) — r/(—l) in the integrand of Eq. (|2.33[) is odd under I — > —I, while 
the other factor, within the curly brackets, is even, the integrand is odd and therefore Eq. 
(12.51) is verified. 
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2.4 Evaluation of the form factor 



Using the trace formula of Eq. (|2.28p , we can write down the individual terms in j„ from 
Eq. (|2,3ip . For this, it is convenient to define the integrals 



I a (pi,P2) 



" Oi 



i (p 2 + 2l- Pl )(p 2 + 2l-p 2 ) 
{(p 2 + 2l- Pl )(p 2 + 2l-p 2 ) 



where 



—S^-rnjMV-vi-l)}. 
These integrals satisfy the relations 



I a (-Pi, -P2) 
I(-Pi,~P2) 



Ia(Pl,P2) , 
~I(Pl,P2) ■ 



(2.35) 
(2.36) 

(2.37) 



(2.38) 
(2.39) 



In addition, since q 2 = 0, both integrals are odd in q if either pi = q or p 2 = q. Using 
these properties, we can write the vertex function in the form 



= 4e/AA' I (4m f —p-p') 



Iu{P,p') + Ia{P,q) ~ Ia(p',q) 



+Pfi 



(p + p , ) a I a (p,p')-q a I a (p,q)-q a I a (p,q)+4m 2 I(p,p') 



(2.40) 



where we have omitted terms proportional to q^ that do not contribute to the amplitude 
due to Eq. (|2.4p . In order to determine the form factor a from this expression, we only 
need to look for the terms proportional to v^, according to Eq. (|2.6|) . Such terms can arise 
only from the integrals 1^, which can be combined in the form 



I/i(p,p') + In(p,q) ~ In(p',q) = -(m 2 - m /2 )4(p,p') , 



where 



IfiiP,^) 



I, 



(2.41) 



(2.42) 



21 ■ q(m 2 + 21 ■ p)(m' 2 + 21 ■ p>) ' 
The result of the integration can be expressed in the form 

I>, p') = Av^ + Bp„ + B'p'^ , (2.43) 
where A, B and B' are Lorentz invariants, and the form factor a is then determined as 



a = 4e\\'(m 2 + m' 2 - 8m 2 )A . 



(2.44) 



The integral defined in Eq. (|2.42p cannot be calculated analytically in the general case, 
as is common for this type of integral that involve the thermal distribution functions. For 
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definiteness here we consider the case of a non-relativistic fermion gas for which a simple 
result can be obtained. Thus, assuming that 



\l \ m f , \l\ <C m f , (2.45) 

then as a first approximation we can neglect I altogether and replace Iq by mj in Eq. 
(|2.43p . As a result, 1^ is proportional to v^, so that the coefficients B and B' are zero, 
while 

A = —r-o WT/ — 7) ^7T / T^W 2 ~ m 2 f )[r](l) - rj(-l)] , (2.46) 

2uj(m? +2m f E)(m' 2 + 2m f E') J (2vr) 3 v fJUKJ n ;J ' v ' 

E' = E-lo. (2.47) 



where 
Therefore, 



m 2 + m' 2 — 8m 2 j 

a = e/AA (n f - n f ) + 2mfE){m , 2 + 2mfE/) , (2-48) 

where nj and are the number densities of the fermions / and their antiparticles, re- 
spectively. 

If the background is charge-symmetric, i.e., contains an equal number of particles and 
antiparticles, the form factor, and thereby the amplitude, is zero as Eq. (I2.48j) clearly 
shows. This result actually follows more generally from the presence of the factor rj(l) — 
r}{—l) in Eq. ()2.37|) . and it holds whether the fermion gas is non-relativistic or not. On 
the other hand, Eq. (|2.48p confirms explicitly that the amplitude is non-zero in general, 
giving a non-zero rate for the radiative decay process. 

Finally, the explicit expression for the differential and total decay rates in this limit 
(non-relativistic gas) are given by substituting Eq. (|2,48p into Eqs. (|2,17p and (|2.20p . and 
remembering the relation given in Eq. (j2.15|) . The final integration that remains to obtain 
the total rate using either Eq. (|2.20p or Eq. (|2.17p . is straightforward but cumbersome. We 
do not pursue this further here since the details are similar to next case that we consider, 
which we treat in full. 

3 Gravitational decay of a neutrino 

3.1 Kinematical considerations 

We now consider the process 

vi(p)^V2(p') + G(q) (3.1) 
where Q denotes the graviton. The amplitude can be written as 

%M = -ine X P* (q)u 2 (p')T XpUl (p) , (3.2) 
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where e xp is the polarization tensor of the graviton, T\ p is the vertex function, and k is 
related to the Newton gravitational constant G by 



(3.3) 



As we have already mentioned, in the calculation of the rate we neglect the effects of 
the medium in the dispersion relations and wave function renormalization of the external 
particles. Therefore, the spinors satisfy the vacuum Dirac equation while the on-shell 
conditions for the graviton polarization tensor are 



eA P (<?) 
ri Xp ex P {q) 

q x t\ P {q) 



o, 
o. 



(3.4) 



The vertex function is constrained by Lorentz invariance and must be transverse to 
the graviton momentum. Subject to these conditions, the general form for the vertex 
function has been obtained in Ref. pp. It was noted there that the vertex function can be 
decomposed into a tensor and a pseudotensor, 



r (T) , r (P) 
1 \ P + 1 \p 75 



(3.5) 



and that for each with I = T, P, the terms that contribute to the vertex function 
with an on-shell graviton are of the form 



r 



CO 



■ v 



+ 



p-q 
q ■ v 



PXPp 



p-q 

1 

q ■ v 



{PV}\P 



-v\v p 



. (P • q) 2 "' ' ' q ■ v ' p-q 
where the a's and the Us are form factors, and 

{AB} Xp = A x B p + A p B x . 



{Pl}\p + {vi}\ P 



(a^+ifV), (3.6) 



(3.7) 



For the present calculation, the on-shell amplitude takes a simpler form. First, the 
transversality of the vertex function implies that we can make the replacement 

eAp -> £A P + X x q p + q x X p (3.8) 

with any X x , which allows us to choose the graviton polarization tensor such that 

0. (3.9) 



v X e Xp 



Second, since the neutrinos are much lighter compared to the charged leptons and the weak 
gauge bosons, we can neglect the masses in calculating the leading term of the amplitude. 
This also enables us to use the massless limit formulas 



U2{p')iui{p) = 0, 
U2{p')il$ui{p) = 0. 



(3.10) 
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Finally, chirality arguments dictate that the terms with or without the 75 in the amplitude 
are not independent. In fact, in the massless limit, the chirality invariance of the standard 
model interactions imply that the terms a 5 and bg cannot be present, while the remaining 
ones actually appear combined in the form 

iM = -iKe*^u 2 (p') [cp^Pvf + dippiv + Pvlti)] Luxip) , (3.11) 

where c and d independent coefficients that must be computed. Eq. (|3.1ip is the most 
general form of the on-shell amplitude in the chiral limit. 

The formulas already given in Eqs. f)2. 12|) and (|2.17p for the angular distribution and 
the decay rate apply in this case as well. In order to evaluate \M\ 2 we use the relation 

u\u x = Li>, (3.12) 

which is satisfied by the neutrino spinor in the massless approximation that we have used, 
and the analogous relation for 112- The graviton polarization tensor can be expressed in 
terms of the spin-1 polarization vector of definite helicity (with s = ±), which are such 
that 

e£ = (0,e s ), (3.13) 

with 

e s -q = 0. (3.14) 
A particularly convenient representation follows by defining the unit vectors 

q = (sin 9 cos (ft, sin 9 sin (ft, cos 9) , 
e\ = (sin (ft, — cos (ft, 0) , 

e*2 = (cos 6 cos (f), cos 9 sin (ft, — sin 9) . (3.15) 
The spin-1 polarization vectors of definite helicity are given by 

e s = — (ei + ise 2 ) , (3.16) 

and the graviton polarization tensor for a definite helicity is 

ef = e^- (3-17) 

In the most general case, the coefficients c and d introduced in Eq. (|3.1ip can have 
different phases, which is a signal of CP-nonconserving effects. However, as is well known, 
such effects do not show up at one- loop when their only source is the mixing matrix. Thus, 
for our present purposes, we can assume that c and d are relatively real. The generalization 
to the the case in which Imc*ci 7^ is straightforward. In this way we then obtain 

\M\> = 2k 2 P 4 sin 4 9 {\c\ 2 E(E - uj) + 4\d\ 2 + 2c* d(2E - u)} (3.18) 

where we have used the relations 

|e s • z\ 2 = sin 2 9 

s 

^|e s -i| 4 = isin 4 #. (3.19) 
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Figure 2: W exchange diagrams for the one-loop contribution to the v e gravitational vertex 
in a background of electrons. The braided line represents the graviton. 

3.2 The diagrams 

We take the medium to be a thermal background of electrons. To one-loop order, T\ p 
is determined by the diagrams shown in Fig. [2j Those diagrams contribute only to the 
gravitational coupling of v e and not to the other flavors, which lead to non-diagonal 
couplings between the neutrino mass eigenstates through the neutrino mixing matrix. 
There is another set of one-loop diagrams involving the Z boson, but they contribute 
equally to the gravitational vertex of all the neutrino flavors v e ,u,,Ti an( i therefore do not 
contribute to the off-diagonal couplings between the mass eigenstates of standard, weak 
SU(2)-doublet neutrinos. They are relevant if the decay process involves a so-called sterile 
neutrino, but we do not consider that possibility here. 

There are additional diagrams in which one or both of the vector boson lines are 
replaced by their corresponding unphysical Higgs partners. As we will see, the leading 
contributions from the W exchange diagrams to the decay amplitude is 0(1/Myy)> and 
therefore the Higgs exchange diagrams are relevant in general. However, for simplicity, we 
will carry out the calculation in the limit m e — > 0, in which case the diagrams involving the 
unphysical Higgs particles do not contribute. Thus, our results are strictly applicable for 
those physical environments in which the electrons in the background can be considered 
to be relativistic. 

There are also diagrams in which the graviton line comes out from one of the external 
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neutrino legs are one-particle reducible and do not contribute to the leading order term 
of the off-diagonal vertex function T\ p . The proper way to take those diagrams into 
account in the calculation of the amplitude for any given process, is by choosing the 
external neutrino spinors Ui(p) to be the (properly normalized) solutions of the effective 
Dirac equation for the propagating neutrino mode in the medium, instead of the spinor 
representing the free-particle solution of the equation in the vacuum. Since the off-diagonal 
vertex function is zero at the tree-level, those corrections yield higher order contributions 
to the decay amplitude, which we neglect. 

Therefore, we will compute the integral expressions defined by the diagrams shown in 
Fig. [21 in a manner that is consistent with the approximations and idealizations that we 
have outlined above. 



3.3 One-loop expressions for the vertex 

We consider situations in which the VF-boson mass is much larger than all the other energy 
scales, so that it is valid to neglect the thermal effects in the W-boson propagator. In 
addition, we can expand the amplitude in inverse powers of the VF-boson mass and keep 
only the leading terms. The one-loop expressions for the neutrino gravitational vertex 
function under the above conditions has been determined in previous works [H [5] , and 
the results that are relevant for our present purposes can be summarized as follows. 
Up to terms of order the vertex function can be written in the form 

T Ap = U el U* e2 [A\ p + G\ p + H Xp ] , (3.20) 

where U is the neutrino mixing matrix. The terms contained in A\ p are the leading order 
terms in the Fermi constant and they were derived in the first two references cited above. 
Those terms do not contribute to the on-shell amplitude, which is most easily seen by 
noticing that those terms are momentum-independent and therefore do not produce any 
contribution of the type shown in Eq. (|3.1ip . 

The remaining terms in Eq. (|3.20[) are the 0(1/Myy) terms, which were determined 
in Ref. [5]. They are momentum-dependent and the relevant ones for the present work. 
The part denoted G\ p contains all the 0{1/My^) terms that are independent of the weak 
gauge parameter £, while H\ p contains the rest, which depend on £. Borrowing the results 
from that reference, we quote below the relevant formulas for the contribution from each 
of the diagrams in Fig. [2] to these two sets of terms. 

Let us consider G\ p first, which contains all the 0(1 /Myy) terms that are £- 
independent, and let us start with diagram (a). The contribution from this diagram 
is given by 



xp 2 J (2vr) 3 

(4) 



m P (i + qMH4)*W Q pV-p') 

(l + q) 2 



(l-q) 2 

1 P L, (3.21) 
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where 

& ( w\v(k) = j^- (k\ v - k^K) , (3.22) 



V Xp (k, k') = - lx {k + k') p + lp {k + k') X 



->v\ P 



-\- jit 



(3.23) 



is the tree-level electron gravitational coupling (in the massless limit), and rj(l) has been 
defined in Eq. (|2.25p . The factor of 5(l 2 ) appears in Eq. (|3.2ip because we have taken the 
limit m e — > 0, as explained earlier. It should be noted that we have retained only the con- 
tribution due to the thermal part of the electron propagator since, as already mentioned, 
the on-shell amplitude vanishes in the vacuum. Similarly, for the other diagrams, 



G 



(6) 



2 



_dH_ 

(2tt) J 



Th-C XpaP (l -p',l-p) + TTaW A SU' - P) + ^"'xp^WoQ ~ P') 
lvl w 



G 



G 



Ap 

(d) 

Ap 



'O'Xpap 



(2tt) 3 
dH 
(2^ 



S^tyfLlhpL^il-p), 



(3.24) 



where the tensors a\ ppu , a' XpfMU , and C\ PPiU that appear in these expressions are related to 
the various gravitational vertices that appear in the diagrams, and are given by 



^Xppv 
a \ppv 

G\ppv (kj k ) 



VXpVpu ~ ^{vxpVpu + Vppmu) , 

VXpVpu ~ (mpVpu + VXuVpp) , 
TJxpiVpvk • k — k p k u ) — Ti pi/ {k\k p + k x k p ) 
+ k u (r]xpk' p + 7] pp k' x ) + k'^rjxukp + r] pv k\) 
- k • k'(r) Xp j} pv + T] Xu r] pp ) . 



(3.25) 



The expression for the Jj^' b > c ' d ^ which contain the ^-dependent terms, are given by 
the same expressions as Eqs. (|3.21|) and (I3.24|) . but with the replacements 



a(4) 

Cx P pu(k, k ) 



A 



(0 

Waf3 



£,k p k v 



£ \jlxpkpk v — kykyXipy — k p k p rj\ v — k\k v r] pp — k p k u rjxp^ • (3.26) 



Before extracting the physical part from these expressions, we consider the conditions 
required by the gauge invariance of the weak and the gravitational interactions. 
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3.4 Weak gauge invar iance 



The gauge invariance of the weak interactions requires that the on-shell amplitude be 
independent of the parameter £. It turns out that the H terms satisfy 



u 2(p')H^Jui(p) = (x = a,b, c, d) , 



(3.27) 



so that the requirement is indeed satisfied. In order to show this, let us consider specifically 
the first term x = a, where 



H 



(a) 

Ap 



' 2 



(0 



Q/-i)V Xp {l,l-q)1/t^ ap {l-p) 



(l-9) s 



+ 



(l + q) 2 



(3.28) 



Using Eq. (|3.26|) . the first term in the square brackets involves the factor f(f—fi)ui(p), 
which reduces to zero when the massless Dirac equation for u% and the delta function 5(l 2 ) 
are taken into account. Similarly, the second term in square brackets contains the factor 
U2{p')(J — $')f , which also reduces to zero for similar reasons. 

Similar arguments hold for the remaining f/-terms, which we do not consider explicitly 
any further. Therefore, the conclusion is that only the G terms, which are independent of 
£, contribute to the physical amplitude. 



3.5 Physical part of the amplitude 

A by-product of Eq. (I3.27P is that the k^k u term in Eq. (13.221) does not contribute to the 
physical amplitude and therefore for the present purposes we can adopt 

^(*) = -3^T^. ( 3 - 29 ) 
w 

in Eqs. (|3.2ip and (|3.24[) . In the evaluation of the expressions the , we can ignore the 
terms that do not contribute to the on-shell amplitude due to the on-shell relations given 
in Eqs. (|3.4p and (|3.10p . Denoting the latter by G"^ we then write 

Gf p = (jOj G'V + , (3.30) 

where the Cr^i contain the terms that survive. A straightforward evaluation of the 
integral expressions then yields 

G Xp^ = i 2J Xpa + P\Jpa + PpJxa] 7° L , 



15 



G 



/(&) 



G 



G 



/(c) 
Ap 

'(d) 



-2 J a Ap - P\Jap ~ PpJaX + PXPpJa 1° L 

+ [A* - PxJa] [(p + - 2p p7 Q ] i 

+ [Jpa " PpJa] [(P + P')"7A " 2 Pxl a ] L 

+ + p% " P • P^] (^Aa7p + ??pa7A)i • 

-p' a (jXalp + Jpalx)L, 



'Ap - ( Aa7p + ^p«7a) ^ • 

where we have introduced the definitions 

dH 



J, 



a ■> 



(3.31) 



(3.32) 



The integral J a \ p cancels when all the terms are added, while the Lorentz structure of 
J a and J Q( g imply that they are of the form 



J ,y 

Jaf3 



A n v a , 

MoVaf3 + A 2 2V a Vf3 ■ 

Thus, substituting Eq. (|3.33p into Eq. (|3.3ip we obtain 

G' Xp = 2A u PxPpi>L + [3^20 - An(p + p') ■ v]{p^ p + p p ~fx)L . 



(3.33) 



(3.34) 



which, together with Eqs. (|3.2U|) and (|3.3U|) . establishes that the amplitude is of the form 
given in Eq. (|3,lip . with 



3.6 Evaluation of the form factors 



U e iU: 2 -V [3^l 2 o - (2E - u)A n ] ■ 



(3.35) 



The definitions in Eq. (|3.33p imply the following expressions for the coefficients An and 

-420, 



A u 
A 20 

which give the integral formulas 

An : 

-420 : 



v x Jx 



\{j\-v x vPJ Xp ) , 



(3.36) 



dH 



(2tt)* 



(/e /e) j 



6 J (2vr) 3 
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d 3 l 

*l\Ue + h) 



(3.37) 



where f e ^ are the distribution functions of electrons and positrons respectively. In the 
massless limit that we have employed these integrals can be evaluated exactly in terms of 
the temperature T and chemical potential fi of the electron gas, to yield 



An 
A 2Q 



— T + — 

12 I 7T 2 



1 

12 



1 



7tt 2 1 
60 + 2^ + 4vr 2 M 



(3.38) 



The explicit formulas for the form factors c and d are obtained by substituting Eq. (I3.38I) 
in Eq. ([3T351) . 



3.7 Decay rate 

Substituting Eq. (|335j) in Eq. (pU8j) . 



|C/el^e2 



,2 / « 



P 4 sin 4 0(C 1 + C7 2 -") 



(3.39) 



where 



Ci = (3A 20 -EA U ) 2 , 

C 2 = a; Aii(3A 20 -EA n ), 



(3.40) 



with wo given in Eq. (|2. 13|) . The decay rate is then obtained by substituting Eq. (|3.39p 
into Eq. (I2.17p . Using Eq. f)2. 15|) to express 9 in terms of u and making the change of 
variable u> = 0J$y, we get 



r 



ojq I g Km 



2ttV \ Myy 



U e iU* 2 \ I, 



where 



l/r 



dy 



Vl-V 2 1 + 



1 2 



(Ci + C 2 y) 



(3.41) 



(3.42) 



The integral / is of course trivial but the resulting formulas are cumbersome. Introducing 



Rr, 



l/r 1 

dy — 

y n 




n+l _ „n— 1 



n = 1 



] n + \. 



(3.43) 



and noting that R- n +2 = Rn, the result can be written in the form 



I = ^2 a nRn , 



(3.44) 



n=l 
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Figure 3: Variation of the rate with the energy of the incoming neutrino for a charge 
symmetric medium. 

where 

ai = -AdVl-V 2 + 2C 2 (3 

a 2 = Ci(7- 3V 2 ) - 8C 2 VT- 

a 3 = -ACWl-V 2 + 2C 2 (1 

a 4 = C x {l-V 2 ). 

It should be noted that, for small V, 

R n = 2V + 0(V 2 ) (3.46) 

and 

Y,a n = 0(V 2 ), (3.47) 

n 

which imply that 

J2a n Rn = 0(V 2 ). (3.48) 

n 

Thus, despite the overall factor of 1/V in Eq. (|3.4ip . the rate vanishes for V = 0, which 
confirms the statements given in the Introduction based on general grounds. 

As we have already mentioned, we have neglected the electron mass in the calculation of 
the amplitude. Therefore, the results that we have obtained are applicable in situations in 
which the electron gas is extremely relativistic. In particular, the following the hierarchical 
relationship 

T, fx > m e > m > m (3.49) 



V z 



V 2 



v z 



(3.45) 
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must hold among the relevant parameters involved, where "T, fj," on the left-hand side 
means "either T or fi, or both". However, the final formulas still depend on some param- 
eters including the initial neutrino velocity and the state of the electron gas. To give an 
idea of how the rate depends on them, we consider two extreme cases regarding the values 
of T and \i in Eq. (|3.38p . For these purposes we write the formula for the decay rate in 
the form 




(3.50) 



with 

1 4 

F(V) = — ^«ni?n. (3.51) 

1 n=0 



3.7.1 Charge-symmetric medium 

We consider first a medium with zero chemical potential. As Eq. (|3.38[) shows, in this case 

A 11 = 0, A 20 = -^T 4 . (3.52) 
Looking at Eq. (|3.4U|) . we find that C2 vanishes in this case, whereas C\ is given by 

C X = i^T 8 . (3.53) 

57600 V ; 

The rate of the gravitational decay is then given by 

'-- m " uj (g 2 Km\U el U; 2 \) 2 f-?L) 8 F u (3.54) 



115200 \ 1 ez, J \M W 

where F% is given by Eq. (|3,5ip but setting C2 = 0. Furthermore, in this case the angular 
distribution of the gravitons has the specific form 

1 dT sin 4 fl 

rdcostf a (i-y C os0) 2 ' 1 j 

which depends on the velocity of the initial neutrino, but is independent of any other 
parameters. 

Fig. [3] shows the plot of F\ as a function of E/m. For low energy neutrinos, F\ — ► 
as E —* m, as already commented. For high energy (E 3> m) neutrinos, 

ai * -<(§)<*, 
a 2 ^ 4Ci , 
a 3 k -4(^]Ci, 



E, 

ot A ~ (— J Ci, (3.56) 
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while 



Ri ~ 21og(-) , 

nn—l / fp\ n ~ 1 

tfn * r - , (n^l) (3.57) 

n — 1 \m/ 



and therefore i*i grows as 

Fl ~ F (ER) = ^ . (3 . 58) 
3m 

The result for the rate remains valid as long as the relationship E <C Mw is maintained. 

3.7.2 Completely degenerate medium 

We now consider the limit T = 0. For this case Eq. (|3.38p gives 

Au = ^, ^0 = "^, (3-59) 



which in turn yield 



C 2 -4o; 



d AE + 3/i 

However, notice that the relations in Eq. (|3.49p imply in this case that 

C 2 



(3.60) 



Ci 



<1, (3.61) 



for all the values of the parameters consistent with the assumptions that we have made. 
Therefore, the decay rate in this case is given by 



where Fx is the same function that appears in Eq. (|3.54p . Therefore in this case the rate 
grows proportional to E for E <C fJL, or E 3 for E 3> fJL. 

4 Concluding remarks 

These two example calculations explicitly confirm the suggestion that angular momentum 
violating processes can occur in a medium which is completely isotropic, even though the 
state of the medium does not change at all in the process. We can understand intuitively 
how this happens as follows. To elaborate this argument we use a radiative decay process, 
like the one described in Sec. [5J as an example. 
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Consider the decay of a particle ip to a particle ip' and a photon in the rest frame of 
the decaying particle. Denoting the total angular momentum of the initial and final states 
by J and J' respectively, in the vacuum we have 

J = L^, + S^, , 

J' = Lip> + L~f + S> + S y , (4.1) 

where L x and S x stand for orbital angular momentum and spin of each particle. Angular 
momentum conservation implies in particular that 

(J- q) = (f ■ q) , (4.2) 

where the brackets in the left and right-hand sides denote the expectation value with 
respect to the initial and final states, respectively. However, since p = 0, the orbital 
angular momentum of the decaying particle is zero and, since and q = —p', neither Lm 
nor L 7 have a component along q. Therefore, conservation of angular momentum implies 

(Sf ■ q) = (V • Q + S y ■ q) , (4.3) 

so that the helicity along the direction of motion of the final particles is conserved in this 
frame. Since the helicity can only be zero for the scalars and ±1 for the photon, the 
process involving scalars is forbidden. 

The diagrams of Fig. Q] in which one of the fermion lines denotes the thermal on-shell 
part of the fermion propagator, physically corresponds to the process 

f(l)+^(p)^f{l)+tl/{j/)+^q), (4.4) 

where / is the fermion in the background medium. Then in place of Eq. (|4.3p . the relevant 
condition is now 

(L f ■ q + S f ■ q + ■ q) = (L f ■ q + S f ■ q + S,^ ■ q + S~f ■ q) , (4.5) 

where the various terms on either side of this equation should of course be added according 
to the angular momentum addition rules. Even assuming that we may drop the spin 
contributions of the background fermions upon averaging over an unpolarized background, 
their orbital angular momentum does not vanish because in the rest frame of the decaying 
particle the fermions appear to be moving in a preferred direction. Therefore, helicity 
is not a good quantum number anymore and the helicity argument does not apply. The 
argument can be extended with minor modifications to the gravitational decay as well. 

It should be emphasized that, while the model of Sec. [2] is a toy model, in Sec. [3] we 
have worked with nothing else than the standard model augmented with the linearized 
gravitational couplings. Thus, the effects calculated in Sec. [3] constitute real predictions 
of the standard model, subject to the approximations we have made. The fact that these 
processes are forbidden in the vacuum by angular momentum conservation arguments, 
implies that, in the medium, their angular distribution and differential decay rates have a 
distinctive form. This could lead to observable consequences in specific physical contexts 
despite the fact that there may exist other competing processes with comparable total 
rates. 
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